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We present a quantitative analysis of the experimental accessibility of the Tonks-Girardeau gas in 
the current day experiments with cigar-trapped alkalis. For this purpose we derive, using a Bethe 
anzats generated local equation of state, a set of hydrostatic equations describing one-dimensional 
5-interacting Bose gases trapped in a harmonic potential. The resulting solutions cover the entire 
range of atomic densities. 
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Physics of trapped atomic gases attracted an enormous 
attention of both experimentahsts and theorists in the 
last decade [0,^. AH atomic systems experimentally ac- 
cessed so far can be well understood within a so called 
mean-field picture where the atom-atom correlations are 
weak and a particular atom can be effectively viewed as 
moving in a mean-field created by others and unaffected 
by the presence of the probe atom. Small deviations from 
the mean-field predictions can be calculated using the 
Bogoliubov approximation |^]. An interesting question 
arrises as what happens beyond the scope of the mean- 
field behavior. For the three-dimensional systems the 
far-beyond-mean-field regime has neither been reached 
experimentally (the JILA experimental group has made 
a significant progress towards this goal 0), nor under- 
stood theoretically. 

For quasi-one-dimensional systems (cigar-shape traps) 
the situation is different. Recall that a one-dimensional 
system of zero-range interacting bosons in a flat-bottom 
box is exactly integrable |5) via Bethe anzats, for all val- 
ues of the coupling strength. The exactly known equation 
of state in a box provides an input for a unified beyond- 
mean-field treatment of trapped quasi-one-dimensional 
atomic gases, and that is what our paper is devoted 
to. The main goal of the paper is to provide quantita- 
tive criteria for an experimental realization of the Tonks- 
Girardeau gas 1^ in atomic experiments, with the main 
emphasis on deviations of the spatial distribution from 
the mean- field Thomas-Fermi profile Q . We paid special 
attention to the intermediate between weak and strong 
interactions regime, as the most realistic from the ex- 
perimental point of view. The extreme case of infinitely 
strongly interacting atoms in a harmonic potential has 
been already investigated in Ref. Also the beyond- 
mean-field effects in the self-correlation function of the 
trapped one-dimensional gas were considered in 

Experimental progress in creating one-dimensional 



atomic gases is also quite fast. One-dimensional guid- 
ing of thermal cold atoms through elongated magnetic 
configurations pO| has already been successfully demon- 
started. A recent experiment shows evidence of trap- 
ping a Bose-Einstein condensate by a shallow magnetic 
gradient in one direction combined with a blue detuned 
hollow laser beam ensuring tight confinement in the ra- 
dial directions. To ensure a truly one-dimensional be- 
havior the atoms must be cold enough to predominantly 
occupy the ground transverse vibrational state. In this 
respect the most promising candidate is atoms in an array 
of ID traps formed in the intersection of two far-detuned 
Gaussian standing waves , obtained in turn by releas- 
ing a full, three-standing-waves optical lattice, where up 
to ~20% ground vibrational mode occupation is already 
experimentally accessible A 2D version of such an 
array (one standing wave) has been already reported in 
literature [ p^ . 

In highly elongated cigar-shaped traps (a;_L ojz), the 
transverse atomic motion is governed by the Hamiltonian 



2m 



(1) 



where m is the mass of the atoms, and uj± and lUz 
are the transverse and longitudinal frequencies of the 
trap respectively. If neither temperature nor interac- 
tion energy-per-particle e (defined below) exceed the 
transverse level spacing huj±, atoms occupy the ground 
mode of this Hamiltonian, and the system becomes ef- 
fectively one-dimcnsional p^ . In particular, the inter- 
particle interaction in the longitudinal direction can be 
well-approximated by an effective two-atom interaction 
potential |l5| 



(2) 



where giD = —h^/flaiY), is an effective one-dimensional 
coupling constant, and oid = (— aj_/2a) [1 — C{a/a±)] 
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is the one-dimensional scattering length defined anal- 
ogously to the three-dimensional case as aiD = 
—d/!^/dkz\k^^Q+, ^{kz) being the scattering phase of 
the even scattered wave. Here and below fl = rn/2 
is the reduced mass, a± = y^h/fluj± is the size of the 
ground state of the transverse Hamiltonian in Eq. (|^), 
a is the three-dimensional scattering length, and C = 

lim.^oo ( Jo ds'/V^ - V\/^) = 1-4603 .... 

Up to terms containing C, the effective one-dimensional 
interaction is a simple projection of the three-dimensional 
zero-range interaction on the transverse ground mode. 
Transverse renormalization effects captured in the C- 
dependent terms become important only for strong con- 
finement, \a\ S> a±. 

In addition to transverse trapping and interparticle in- 
teractions, the atoms are also subject to weak residual 
longitudinal trapping, which is usually represented by a 
one-dimensional harmonic potential of a frequency w^. 
The effective one-dimensional Hamiltonian for N trapped 
atoms is thus 
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where 
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is the well-known Hamiltonian for a onc-dimcnsional 6- 
interacting Bose gas in a flat-bottom box. This Hamilto- 
nian can be diagonalized via Bethe anzats j^, and thus 
the equation of state of such a gas is known exactly for 
all densities and temperatures . At zero temperature 
the energy per particle e(n) is given through 



(5) 



where the dimensionless parameter 7 — 2/n\aiD\ is in- 
versely proportional to the one-dimensional gas param- 
eter n IoidI, n is the one-dimensional number density of 
particles, and the function 6(7) is given by 



e(7) 



^y^^.g(x|7)x dx, 



(6) 



Functions g {x\^) and A(7) are solutions of the Lieb- 
Liniger system of equations 
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A(7) =7 j ^ gi^l) dx . (8) 

Limits of small and large gas parameter n |aiD| can be 
expressed in closed form (see S): 



as n IoidI — > 0, e{n) 



6m 
1 



(9) 



as n IoidI — > 00, e(n) ^ - 51D ?^ ■ (10) 

The low-density limit (^ corresponds to the case of in- 
finitely strong interactions; the corresponding atomic sys- 
tem is usually referred as a gas of impenetrable bosons or 
Tonks-Girardeau gas Notice that the expression for 
the energy-per-particle formally coincides with the one 
for free fermions: this is a manifestation of the Fermi- 
Bose duality in one-dimensional systems . The oppo- 



site limit of high densities (10) represents (if multiplied 
by the number of particles and integrated over volume) 
the Thomas-Fermi energy functional first introduced by 
Bogoliubov 1^. It is usually viewed as the thermody- 
namic limit of the Gross-Pitaevskii energy functional 
p8| , derived in turn on the basis of the mean-field ap- 
proximation. It should be mentioned that in the three- 
dimensional case the mean-field approximation is valid 
for low densities contrary to the one-dimensional case 
where it is valid at high densities instead. 

We now introduce the classical hydrodynamical ap- 
proximation: we suppose that the trapped gas at each 
position z is in local thermal equilibrium, with local en- 
ergy per particle given through Eqn. (|^). At zero tem- 
perature the hydrodynamic equations of motion read 
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—V + v—v = -(1/m) — (0(n) + Viz)) 
ot oz oz 



(11) 



where v is the atomic velocity, V{z) = 
potential energy per particle, and 

T— d 
(j){n) = (1 + n— )e(n) 



moj^z^ /2 is the 



(12) 



is the Gibbs free energy per particle obtained using the 
Eqns. II and |lll-|| and tabulated in Ref. |l|l. The va- 
lidity of the hydrodynamical approximation requires the 
typical energy per particle e ~ mw^/2 -I- -I- ^ be much 
higher than the longitudinal quantum level spacing TilOz- 
In what follows we will be interested in a steady-state 
solution of the system (pT]). The first integral of the sta- 
tionary version of this system obviously reads 



0(n) + V{z) = ^ 
n = 



for \z\ < R 
for \z\ > R 



(13) 



where /i can be proven to be the chemical potential of 
the system, fixed by a normalization condition 



n[z) dz — N 



(14) 



Here the atomic cloud radius R is given by V{R) — ^. 
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FIG. 1. Discrepancy between the cloud size computed from 
the exact equation of state and from the Thomas-Fermi pre- 
diction, as a function of the (single) governing parameter rj 
( p3| ) shown in the upper a;-axis. As an example of a "physical" 
variation of the parameter ri we consider A'' = 200 confined in 
a cigar-shape trap of transverse frequency i'± = 180 kHz ^'^''Cs 
atoms, whose scattering length is magnetically tuned [20] to 
a — -1-120 aBohr- Variation of the parameter rj is supposed 
to be realized through the variation of the longitudinal fre- 
quency (lower a:-axis). For the three representative points 
marked in this Figure, we computed the density profiles from 
the exact, Thomas-Fermi and Tonks-Girardeau equations of 
state. Results are shown in Fig. ^. 

In the limit of low density (||) the density profile is a 
square root of parabola B 



n{z) = nJ^„„ks 



1 - 



^Tonks 
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(15) 



for z G [— i?Tonks, -Rxonks], with n = elsewhere, where 

n'Ton^,^[i2/n^)Nimu;Jhf^' ; (16) 
RToni.s = [2N{h/Lo,m)f^ . (17) 

The opposite limit ( p^ reproduces the familiar Thomas- 
Fermi parabola 0|. In this limit, it is easy to show that 



n{z) = n^p ^ 



1 - 



p2 



for z e [— i?TF, ^tf], with n = elsewhere, where 



i?TF 
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(18) 

(19) 
(20) 



In general, the rescaled density s — n/n^p satisfies the 
rescaled version of the system (O) and (111), 
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rj s{x) 

b 



1 2 



s{x) dx — — 
o 



(21) 
(22) 



with X = z/Rtf is the dimensionless coordinate, and 
fil) = [3e(7) — 7(ie/d7] is the dimensionless Gibbs 
energy (see |l^). We see that this system has only one 
governing parameter, 



77 = n^p IaidI 



(9/2)1/3 faimNLU.il -Ca/aj_)^^'^^^ 



(23) 



whose physical meaning is the Thomas-Fcrmi estimate 
of the one-dimensional gas parameter in the center of 
the trap. It is easy to show that the limit of small 77 
corresponds to the Tonks-Girardeau limit, and the limit 
of large 77, to the Thomas-Fermi limit: 



77 <C 1 — » Tonks-Girardeau profile (15) 
77 ^ 1 — > Thomas-Fermi profile (nSl) 



(24) 
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FIG. 2. Comparison of density profiles produced by 
Thomas-Fermi, exact, and Tonks-Girardeau local equations 
of state for the points indicated in Fig. |l|. (a) rj = 0.07, 
= 10 Hz; Tonks-Girardeau regime, (b) 77 = .4, fz ~ 130 Hz; 
intermediate regime, (c) rj = 9, — 13 kHz; Thomas-Fermi 
regime. The number of atoms, three-dimensional scattering 
length, and transverse frequency are the same as in Fig. ^. 

The exactly known equation of state Eqs. (H) allows 
one to precisely map the transition between these two 
limiting behaviors. A convenient way to detect, both 
experimentally and numerically, the degree to which the 
system is in one or the other limiting regime is to compare 
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the size of the cloud of trapped atoms, R computed from 
the Eqs. ( p^ , p^ , to the Thomas-Fermi prediction i?TF 
in Eq. (20). This is shown in Fig. ^ Density profiles for 
three representative values of the governing parameter 
77 are plotted in Fig. |[ As the Figs.|l|-|| make clear, the 
transition between the two limiting regimes happens over 
the range of the governing parameter 77 from about 0.1 
to 10, reaching 75% of the cloud size discrepancy for the 
former. 

Let us summarize the requirements for a successful ob- 
servation of the zero-temperature Tonks-Girardeau den- 
sity profile (p^): 

(a) The couphng constant giD (see (§)) must be posi- 
tive. This condition is equivalent to < a < a±/C; 

(b) To ensure the one-dimensional behavior of the sys- 
tem the energy per particle must not exceed the 
transverse level spacing: in the case of the Tonks- 
Girardeau gas this condition leads to huj^N <C 

(c) The number of particles must be well below 
the Tonks-Girardeau vs. Thomas-Fermi boundary 
which requires 77 <C 1, where rj is given by (p3|); 

(d) It turns out that not for every set of trap and atom 
parameters the energy per particle at the Tonks- 
Girardeau vs. Thomas-Fermi boarder (77 ~ 1) is 
consistent with the hydrodynamical approxima- 
tion, i.e. exceeds the longitudinal level spacing. To 
avoid that one must require /m |aiD|^ ^ Tti^z or 
a ;» {a\/az){'^ — Ca/a±), where ~ (?i/77i(jjj_)^^^ 
is the size of the longitudinal ground state. If the 
above condition is violated, then, as the number 
of particles decreases, the system will pass from 
the Thomas-Fermi regime to the ideal gas directly, 
skipping the Tonks-Girardeau regime; 

(f) One may require also, that when already in the 
Tonks-Girardeau regime (77 <C 1), the number of 
atoms must nevertheless be high enough to ensure 
the hydrodynamical, non-ideal gas behavior. It 
turns out that this requirement leads to a trivial 
condition of having more than one atom in the sys- 
tem: > 1. 
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